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Synthetic nanoscale complexes capable of mechanical movement are often studied theoretically 
using discrete-state models that involve instantaneous transitions between metastable states. A 
number of general results have been derived within this framework, including a "no-pumping theo- 
rem" that restricts the possibility of generating directed motion by the periodic variation of external 
parameters. Motivated by recent experiments using time-resolved vibrational spectroscopy [Panman 
et al., Science 328, 1255 (2010)], we introduce a more detailed and realistic class of models in which 
transitions between metastable states occur by finite-time, diffusive processes rather than sudden 
jumps. We show that the no-pumping theorem remains valid within this framework. 



I. INTRODUCTION 

Synthetic molecular complexes capable of mechanical 
movement, often simply called molecular machines, have 
caught the attention of researchers in the past decade for 
their potential applications in nanomechanical devices 
[TJ [2]. Numerous designs have been proposed and ex- 
perimentally realized, including unidirectional rotation 
of molecular components [3] [4], molecular translation 
on tracks [5] and surfaces [6], a programmable DNA 
nanoscale assembly line [7] and a single molecule electric 
motor [8 j. The operating principles of these molecular 
machines differ from those of their macroscopic counter- 
parts: inert ial effects are negligible, whereas thermal fluc- 
tuations play a dominant role, making their description 
inherently stochastic. 

A number of recent theoretical studies have focused on 
the response of molecular machines to the time-periodic 
variation of external parameters, a mode of operation 
often called stochastic pumping. In this setting, the cur- 
rents induced in the molecular machine have been shown 
to have geometric contributions [9HT3] analogous to the 
Berry and Aharonov-Anandan phases in quantum me- 
chanics [14[|T5]. Moreover, in Ref. [16 an exact expres- 
sion for the average pumped current was derived, and 
was used to establish a no-pumping theorem (NPT), a 
restrictive condition for generating directed motion in 
stochastic pumps. These results were further studied and 
extended in Refs. [T7H24] . Most of this theoretical work 
has involved discrete-state models in which the molecular 
machine makes instantaneous, thermally activated tran- 
sitions among a set of metastable states. Such models 
are conveniently visualized using a network representa- 
tion (see e.g. Fig. |4|, where the nodes of a graph depict 
the metastable states and the edges indicate the allowed 
sudden transitions [25j [26] . 

In reality, a transition between two metastable states 
of a molecular machine involves mechanical motion, and 
therefore can not be instantaneous. Recent experiments 
[271 128] using time-resolved vibrational spectroscopy to 
study the movement of a molecular machine between 
two docking stations provide evidence that this motion 



is described more accurately as a rapid, one-dimensional 
random walk than as an instantaneous jump. This ran- 
dom walk description can capture the effects of physical 
separation between the docking stations. Motivated by 
this observation, in this paper we introduce a model of 
molecular machines in which the system makes diffusive 
(rather than sudden) transitions. The metastable states 
continue to be depicted by the nodes of a graph, but dur- 
ing a transition from one state to another the machine 
evolves diffusively along the connecting edge. Thus the 
machine now lives in a hybrid state space, consisting of 
both discrete and continuous components, represented by 
the nodes and edges of the graph. 



Because our hybrid model aims to provide a more de- 
tailed microscopic description of the dynamics of molec- 
ular machines, it is natural to wonder whether the theo- 
retical results mentioned earlier, which were derived un- 
der the assumption of instantaneous state-to-state transi- 
tions, remain valid within the framework we now propose. 
In this paper we will focus specifically on the no-pumping 
theorem (NPT) of Ref. [16] , which was observed in actual 
experiments [29], in fact prior to its general theoretical 
formulation. We will show that the NPT remains exactly 
valid within our hybrid model. 



The no-pumping theorem was experimentally observed 
in Ref. [29] for a [2]catenane complex, a rotatory molec- 
ular machine. The [2]catenane also provides a minimal 
setup for illustrating key elements of our hybrid frame- 
work and its analysis. The first part of the text is there- 
fore devoted to the discussion of this particular system. 
We then generalize our framework by extending it to 
more complex molecular machines. The organization of 
the paper is as follows: Sec. [TT] i ntroduces the hybrid 
model for the [2]catenane; Sec. Ill imposes detailed bal- 
ance; Sec. |IV| provides the proof of the no-pumping theo- 
rem; and Sec. [V| generalizes these discussions to arbitrary 
systems. 



2 




(a) 



1 2 3 



(b) 



FIG. 1. A [2]catenane and its representations, (a) Block 
diagram, (b) Hybrid model. Dashed line indicates periodic 
boundary condition. 



II. HYBRID MODEL OF A [2] C ATEN ANE 

A [2]catenane is a supramolecular complex consist- 
ing of two mechanically interlocked rings, called macro- 
cycles. In the experiments of Ref. [29 these macro- 
cycles were unequal in size, and the smaller macrocy- 
cle had three binding sites, or stations, on the larger 
macrocycle. In Fig. [TJa) these stations are represented 
schematically by numbered boxes. This system has three 
met ast able states, corresponding to the presence of the 
smaller macrocycle at station 1, 2 or 3. 

Consider now the thermally activated motion of the 
small macrocycle along the large macrocycle. In the dis- 
crete state model the state space of the system is depicted 
by three points, among which the macrocycle makes in- 
stantaneous transitions, representing sudden jumps from 
one station to another. In our hybrid model, shown in 
Fig. [ljb) , the discrete states are connected by continuous 
line segments, or tracks, along which the system performs 
diffusive motion during each transition. The dashed line 
indicates periodic boundary conditions. 

Let Pi(t) denote the probability to find the system in 
station z at time t, and let Pi(x,t) be the probability 
density to find the system at a position x along track i 
at time t. In our notation, a given track is designated by 
the same index as the station on its left; x specifies the 
distance along a track; and for simplicity we assume each 
track to be of length I. See Fig. [2|a) for an illustration. 
Because of the periodic nature of the state space we make 
the identifications: z + l = lifz = 3 and i — 1 = 3 if 
i = l. The total probability is normalized to unity: 



Y^Pi{t) + Y. dx Pi (x,t) = l. 

-_i JO 



(1) 



Here, ns = 3 is the number of stations and nr = 3 is the 
number of tracks, but in general these need not be equal 
(see Sec.0. 

We now specify the dynamics of our model. When the 
system is in station z, it has a probability per unit time 
cti to make a leftward transition to the location x = I on 
track i — 1, and similarly a probability rate to make 
a rightward transition to the location x = on track z; 
see Fig. j2^b) . When the system is on one of the tracks, 
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FIG. 2. Elements of the hybrid model, (a) Probability at 
station i is denoted by Pi (t) , and probability density on track i 
by pi(x, £), where distance x is measured from station i. Each 
track is assumed to be of length /. (b) Rate parameters of the 
model: a* and f3i denote the transition rates from station i 
to tracks (i — 1) and z, respectively; Iji and ISi are associated 
with the reverse transitions, as described in the text. 



it performs diffusive motion with a fixed diffusion con- 
stant D, with reflective (hard- wall) boundary conditions 
at x = and x = I. Upon reaching either end of the track 
the system might jump into the adjacent station. These 
jumps are characterized by probability rates ji (for tran- 
sitions from track i — 1 to station z) and S{ (from track i to 
station z). More precisely, the probability per unit time 
for the system to make a transition from track i — 1 to 
station z is given by the product ljiPi-i(l,t), where the 
factor I is introduced on dimensional grounds, so that the 
parameter 7^ has units of a probability rate i.e. (time) -1 . 
A similar expression holds for transitions from track i to 
station z; see Fig. [5|b). 

These transitions give rise to a flow of probability be- 
tween stations and adjacent tracks. The net current from 
track i — 1 to station i is given by 



J-+i(t) = ljiPi-l(l,t) -OtiPi(t), 

and that from station i to track i by 

J ir +(t) = p i P i (t)-l5 iPi (0,t). 



(2) 



(3) 



(The subscripts "— z" and "z -V indicate rightward 
probability current into and out of station z, respec- 
tively.) The diffusive current at position x along track 
z is 



J?(x,t) = -D-^ Pi (x,t), 



(4) 



and the reflective boundary conditions imply that 
Jf(0,t) = jf(l,t) = 0. 

These currents generally lead to changes in the proba- 
bility distribution. The rate of change of the probability 
to find the system at station i is the difference between 
the incoming and the outgoing currents, 



dPjjt) 
dt 



(5) 
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and that of probability density along track i obeys a dif- 
fusion equation with a source and sink: 



dpi(x,t) 
dt 



d_ 
dx 
-6(x 



Jf(x,t) 
•Z)J_> 



+ 5(x- 



O)J^(i) 



(6) 



Eqs. [5] and [6] form a set of six coupled, linear equations 
(taking i = 1, 2, 3) which collectively constitute the mas- 
ter equation describing the stochastic evolution of the 
system. 



III. CONSTRAINTS IMPOSED BY DETAILED 
BALANCE 
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FIG. 3. Free energies and barriers. Station i has energy Ei, 
and track i energy e^. Track (i — 1) and station i are separated 
by barrier Bi^. Station i and track i by Bi, r . 



Since our model is meant to represent a system im- 
mersed in a thermal reservoir, the dynamics described 
by our master equation should have the property that 
when the rate parameters c^, ft, ji and Si are held fixed, 
the system relaxes to a state of equilibrium in which all 
currents are zero. This condition of detailed balance im- 
poses constraints on the rate parameters. We now ex- 
plore these constraints, and show that they allow us to 
rewrite the rate parameters in the suggestive form given 



by Eq. 11 below. 

Let Pf q and p e i q {x) denote, respectively, the station 
probabilities and track probability densities in the equi- 
librium state. According to the condition of detailed bal- 
ance, the currents appearing on the left sides of Eqs.[2|[3] 
and 2] vanish when these values are substituted into their 
right sides. This leads to the relations, 



I6ip?(0) 
d_ 
dx 



0. 



: T]exp[-B ijL } 
r]exp[-B i>R ] 



(7) 
(8) 

(9) 



Eqs. [7] and [8j together with an arbitrary frequency scale 
77, define the dimensionless parameters B^l and £^ ?j r, 
while Eq. [9] implies that the equilibrium probability den- 
sity is uniform along each track. Introducing the dimen- 
sionless parameters 



Ei = - In P? q and e* = - In (Zpf ) 
now allows us to rewrite Eqs. [7] and [8] as follows: 

Oi = rie- Bi ' L+Ei , Pi = rie- Bi > R+Ei , 



(10) 



(11) 



Because these expressions resemble the Arrhenius form 
for chemical reaction rates, it is natural to interpret Ei 
(or €i) as the energy of the small macrocycle when it is at 
station i (or on track z); and Bi^ (or B^r) as the height 
of the barrier that separates station i from the track im- 
mediately to its left (or right). These energies are given 
in units of kT, where k is Boltzmann's constant and T the 
absolute temperature of the environment. We note that 



this interpretation, illustrated in Fig.[3j is convenient but 
not strictly necessary: the NPT that we derive below is 
an exact result stated in terms of the well-defined quan- 
tities Ei, ei, B iy L and B iy R, regardless of whether or not 
we choose to interpret these as energies and barriers. 



IV. STATEMENT AND PROOF OF NPT 

Stochastic pumping corresponds to the periodic varia- 
tion of the transition rates {^,^,7^,^}, subject to the 
constraints imposed by detailed balance. In the energetic 
picture introduced above, this translates to the periodic 
variation of the energies and barriers {Ei, e^, B^l, Bi f>}, 
that is, Ei(t + t) = Ei(t), etc., where r is the period of 
the pumping. Under these conditions the system relaxes 
to a unique periodic steady state, 

Pf\t + r)=Pr{t) , p p i s (x,t + r)=p p i s (x,t), (12) 

characterized by time-periodic currents passing through 
the stations, J^(t) and Jf_^(t), and along the tracks, 
Jf' ps {x,t). We are interested in the integrated currents, 

^{x)^ J^jf^{x,t), (13) 

where J r denotes an integral over one period of pump- 
ing. Here, represents the net flow of probability 
from track i — 1 into station i over one pumping cycle, 
and <J>^ and &f ps (x) have similar interpretations; the 
integrated currents thus measure the extent to which the 
pumping of the energies and barriers drives a non-zero 
current around the cycle depicted in Fig. [TJb). Physi- 
cally, these currents measure our ability to generate di- 
rected mechanical motion of the small macrocycle around 
the large macrocycle, by the periodic variation of external 
parameters, with positive and negative <£'s corresponding 
to clockwise and counterclockwise motion, respectively; 
see Fig. [TJa). 
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These considerations apply to the time-periodic 
pumping of any combination of the parameters 
{Ei,€i, Bi L, Bi f>}. In the subsequent analysis, however, 
we will assume that the track energies remain constant 
with time, while the station and/or barrier energies (the 
E"s and B's) are varied periodically. Thus we treat the 
tracks as fixed conduits for diffusive motion from one 
station to another; this is in keeping with the relevant 
experimental studies [27j [29] , which did not include any 
time-dependent track energies. The no-pumping theorem 
(NPT) that we now prove states that, in order to gen- 
erate non-zero integrated currents, we must vary some 
combination that includes both station energies (the E^s) 
and barrier energies (the -B^l's and/or B^rs). In other 
words: (1) if we vary only the barrier energies, while 
keeping the station energies fixed, or (2) if we vary only 
the station energies, while keeping the barrier energies 
fixed, then in either case all the integrated currents will 
be zero. To prove the NPT we now consider these cases 
separately. 

The first case is easy to analyze, and intuitively plausi- 
ble. Let the term instantaneous equilibrium distribution 
denote the equilibrium distribution corresponding to the 
current parameter values, at any instant in time. By 
Eq. [lOj this distribution depends only on the state en- 
ergies and not on the barrier energies. Thus when the 
E^s and e^'s are held fixed, the instantaneous equilibrium 
distribution {P^ q ^p^ q (x)} also remains fixed. Moreover, 
by Eqs. [5] and [6] this time-independent distribution is a 
stationary solution of the dynamics, even if the barrier 
energies are varied with time. At the same time, any 
time-periodic pumping protocol uniquely determines a 
time-periodic solution of the dynamics. It follows that 
when the state energies are held fixed and the barrier 
energies are varied periodically, then the unique periodic 
steady state of the system is simply the time-independent 
instantaneous equilibrium distribution. In this state, all 
the instantaneous currents are zero, and therefore the 
integrated currents also vanish. 

Now consider the situation in which the station en- 
ergies are varied periodically in time, Ei(t + r) = Ei(t), 
and the barriers {B iy L, Bi^} (together with the track en- 
ergies €i) are kept fixed. The NPT then follows from a 
combination of two conditions: the detailed balance con- 
straints, Eq. [TT] and the periodicity of the probability 
distribution, Eq. [l2j as we now show. 

Let us first explore the consequences of the detailed 



balance constraints. Combining Eq. 11 with the expres- 
sions for the instantaneous currents, Eqs. [2] and [3J we 
derive for each station z, 



periodic steady state. Summing both sides over z, we get 



The definition of the diffusive current, Eq. |4j implies 



(15) 



D 



&xjf(x,t)= Pi(l,t)-pi(p,t) 



which combines with Eq. 15 to give 



e B " L J^i(t)+e B ^Ji- 



>(t) + -e 



Jo 



dxJ?(x,t) 



(16) 



= 0. 



(17) 



(See Refs. [2TJ [30] for the analogue of this relation in 
the context of a discrete state model.) If we now assume 
the system has reached a periodic steady state, and we 
integrate this relation over one period, we get 



£ 



ps 



rjl 



Jo 



dx$f ps (x) 



0. 



(18) 



Now we explore the implications of the periodicity of 

Since the probabil- 



the probability distribution, Eq. 12 



ity to find the system in station i returns to the same 
value after each period, the integrated current that en- 
ters that station from the left must be balanced by the 
integrated current that exits from the right: = 3>f^. 
This value is in turn equal to the integrated current enter- 
ing track i from the left. Along track i the integrated cur- 
rent &f ps (x) must be the same for any two points x\ and 
X2, otherwise there would be a net accumulation or deple- 
tion of probability in the interval between those points, 
over each period. Proceeding in this manner around the 
entire circuit we conclude that the integrated current is 
uniform all along Q 



(19) 



It follows immediately from Eqs. [18] and [19] that all in- 
tegrated currents are zero. Thus the NPT is established 
for the three-state model depicted in Fig. [I] 



V. GENERALIZATIONS 

We now generalize our discussion along two different 
directions. First, following Ref. [27] . we allow for a spa- 
tially non-uniform (but still time-independent) free en- 
ergy landscape along each of the tracks. Secondly, we 



e Bi - L J^i(t)+e Bi > R Ji-+(t) = ril[e ei - 1 p i - 1 (l,t)-e*p i (p,t)]. 

(14) 

Note that the superscript ps does not appear here, as this 
relation is valid whether or not the system has reached a 



1 A derivation of Eq. |l9| directly from the master equations requires 
a more careful treatment, in which the source and sink terms in 
Eq. [6]are displaced slightly from the track-ends. As this does not 
contribute conceptually to the main line of the proof we present 
the analysis in Appendix [A] 



move beyond the simple three-station, three-track net- 
work shown in Fig. [I] and extend our model to encom- 
pass an arbitrary, finite network of stations and tracks. 
The analysis involved in these more general situations is 
similar to that presented in Sec. [TV] therefore to avoid 
repetition we sketch only the key ideas in the following 
discussion. 

First we allow a nonuniform energy landscape Vi(x) 
along each track z, instead of constant e*, again in units 
of kT; this leads to the expressions, 



J?(M) 



-D 



dpi(x,t) dVi(x) 



dx 



dx 



Pi(x,t) 



(20) 



-Bi.L+Vi-iCO 



ji = rje^ ,r *-^ , 5 i = 'qe-" i > R+Vi WK (21) 

(compare Eqs.|4]and 11 respectively). When all the state 
energies are held fixed and only the barrier energies are 
varied with time, the arguments presented earlier apply 
here without modification, and we can conclude that all 
currents vanish in the steady state. When instead the 
barrier energies are fixed and the station energies are 
varied periodically, in place of Eq. dl4|) we have 



-B t , R +Vi(0) 



e Bi > L J^(t) + e Bi > R j£(t) 



= rjl 



oVi- 



i(0 p£i(M) 



(22) 



which leads to a modified version of Eq. (18), 



'•**&+7^ [ l dxe Vi W$f ps (x) \ =0. 
D Jo . 



(23) 

The periodicity of the probability distribution, Eq. [l2j 
again implies a uniform integrated current, Eq. [l9j The 
combination of Eqs. [19] and [23] in turn immediately im- 
plies that all integrated currents vanish, and thus the 
NPT is established. 

Now consider a more general network, composed of 
ns stations and tit tracks,, which need not be equal. 
We assume this network is connected: for any pair of 
stations i, j, there exists at least one path - a sequence 
of alternating tracks and stations - that connects station 
i to station j. (This assumption is not restrictive, as 
any finite network can be decomposed into two or more 
connected networks, each of which satisfies the NPT.) 
Fig. (4^a) illustrates a connected network, with ns = 4 
and ut = 5. The [2]catenane model analyzed in Sections 
[TT] - [IV] had only one cycle - a closed loop of stations 
and tracks (see Fig. [I]) - but in the more general case 
considered here there can be more than one cycle, as 
shown in Fig. [4^a). 

We will now use Roman indices i G {1, ... 715} to de- 
note stations and Greek indices a G {1, . . . n^} to denote 
tracks. To each track a we assign an arbitrary direc- 
tion, such that the station at one end of the track is 
viewed as the origin o(a) and the station at the other 
end is the destination d(a), as illustrated in Fig. [4^c). 
These assignments are made independently for the vari- 
ous tracks, and a given station can simultaneously serve 
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(a) 



Energy 



d o(x) 




o(a) or : d(a) 



(c) 



FIG. 4. Aspects of generic networks, (a) A network with 
ns = 4 stations and nr — 5 tracks, (b) Notations for a generic 
network, as described in the text. Although we show a station 
connected to two tracks, in general a station can be adjacent 
to as few as one and as many as nr tracks, (c) A single cycle 
within a network. Arrows indicate the (arbitrarily assigned) 
directions of the tracks, and the origin and destination of track 
a are explicitly labelled. 



as the origin of some tracks and the destination of others. 
Station and track energies are denoted by Ei and V a (x), 
respectively; the barrier between station i and track a 
is given by Bi iOL \ the integrated current from track a (or 
station i) to station i (or track a) is given by 3>^L^ (or 
^ff> a ); and the integrated diffusive current on track a 
along the preassigned direction by ^ ps (x). These con- 
ventions are illustrated in Fig. [4^b), where superscripts 
"ps" have been omitted to avoid clutter. 

As earlier, under the time-periodic variation of station 
and/or barrier energies the system settles into a unique 
periodic steady state. Moreover, when all the state en- 
ergies are held fixed and only the barrier energies are 
varied with time, then the entire system simply relaxes 
to a state of thermal equilibrium, in which all currents 



6 



vanish. To establish the NPT, it remains to show that 
the integrated current vanishes when the barrier energies 
are kept fixed and the station energies are varied period- 
ically. 

The periodicity of the probability distribution implies, 
for each station i, 



{a\d(a)~- 



,,ps 



■i} 



{a\o(a) = 



d,ps 
a' 



(24) 



i} 



the sum on the left represents the net integrated current 
into station i from all tracks a for which it is the des- 
tination, and the sum on the right is the net integrated 
current out of the station into all tracks to which it is the 
origin. For each track a we have 



o(a)- 



^ ps (x) = & 



>d(a) 



(25) 



(compare with Eq. 19), which ensures that probability 
neither accumulates nor depletes anywhere on the track, 
with each cycle. Eqs. [24] and [25] have an interesting con- 
sequence: if there exists a non-zero integrated current in 
the system, then it must be part of a cycle (of alternating 
stations and currents) along which all of the integrated 
currents point in the same direction. The intuition is 
straightforward: to prevent the systematic accumulation 
of probability within the network, current must flow in 
a circle. We now formalize and establish this statement, 
and then use it to prove that all integrated currents must 
be zero (when the barriers are fixed and the station ener- 
gies are varied periodically), following arguments similar 
to those presented in Ref. [2T] , 

Without loss of generality, suppose that a particular 
track a supports a positive integrated current: $ J ps > 0. 
(Equivalent arguments would apply if the sign were neg- 
ative.) Eq. [25] then implies positive integrated currents 
^o(a)^a an d ^a%d(a)* ^ s P r °bability cannot deplete on 
station o(a) over a complete cycle, o(a) must have neigh- 



boring track(s) (3 such that ^% Q ^ > 0. Similarly, to 
avoid accumulation of probability, d(a) must have neigh- 
boring track(s) 7 such that 3>d( a )_). 7 > 0- The periodic 
conservation of probability in turn establishes the direc- 
tionality of the integrated currents along these tracks: on 
track /?, the integrated current must flow toward station 
o(a), and on track 7 it must flow away from station d(a). 
Continuing in this manner, we now construct a set D(a) 
of stations and tracks to which there is a positive flow of 
current from a; this will consist of d(a),7, so on. Simi- 
larly, we construct set S(a) of stations and tracks from 
which there is a positive flow to a; this will consist of 
o(a),/3, so on. In order to prevent the accumulation of 
probability in the former set and its depletion in the lat- 
ter with each complete cycle, D(a) and S(a) must have 
a common element. This implies the existence of a cycle 

c = a -» d(a) -» 7 -» . . . -)> ft -» o(a) — )> a, (26) 

along which all integrated currents flow in the same di- 
rection. 



Now recall that each track in our network has been as- 
signed a direction, pointing from its origin to its destina- 
tion. By assumption, for track a this direction is parallel 
to the direction of probability flow around the cycle c, 
indicated by Eq. [26] However, since the assignment of 
track directions is arbitrary, each of the remaining tracks 
in the cycle (/?,-•• 7) might be directed either parallel or 
anti-parallel to the flow along the cycle. Let us therefore 
introduce a factor = ±1, defined for every track /i in 
the cycle c, such that = +1 (or —1) if track /i is ori- 
ented parallel (or anti-parallel) to the flow in the cycle. 
Then ^ d ^ ps = s^ d ^ ps > for each track fi in the cycle 
c, and V ps f , = s„<$> ps f , and , ee s^* , 

are positive as well, by Eq. [25] The analogue of Eqs. [18] 



and 23 for the cycle c is given by 



E 



+ t! / ^e v ^ x) ^ ps (x) 
D Jo 



0, (27) 



where the sum is taken over all tracks \i in cycle c. Since 



all the ^'s are positive, Eq. 27 cannot be satisfied. Hence 
our starting assumption, the existence of a positive in- 
tegrated current along track a, must be invalid. This 
establishes the no-pumping theorem. 



VI. CONCLUSION 

Motivated by the experiments of Refs. [2TJ [28] we have 
introduced a hybrid class of models for artificial molec- 
ular machines, combining elements of both discrete and 
continuous models. Within this picture, we have estab- 
lished the validity of the no-pumping theorem [161 [17], 
which had earlier been derived for entirely discrete and 
entirely diffusive models. We illustrated our deriva- 
tion in detail for the [2]-catenane complex studied in 
Ref [29], and then generalized it to more complicated 
situations. It remains to be investigated whether other 
results [H [18J [201 H2 [30] that have been derived for 
purely discrete or purely continuous models, apply to our 
hybrid model. 
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FIG. 5. Source and sink terms are moved away from the 
track-ends by a parameter < e < 1/2. 



Appendix A: Derivation of Eq. |19| 

Here we establish Eq. [19] directly from the master 
EqsJE] and |6j evaluated in the periodic steady state, 
Eq. [12] We begin by assuming that the source and sink 



for each station i. Equivalently, J dt dPf s (t) /dt = 0, 
which from Eqs. [5] and [13] implies 

® P 4i = ( A2 ) 

This is the first part of Eq. [l9j Similarly, for each track 
i we have J dtdp p i s (x 1 t)/dt = 0, which combines with 
Eq. [AT] and jl3] to give 

^t ps (x) = S(x - e)*£ -5(x-l + eW4 i+1 . (A3) 

From the reflective boundary conditions, jf(0,t) = 
Jf{l,t) = 0, we have 



*?* fl (0) = $?' pfl (Z)=0. 



d,ps / 



(A4) 



terms in Eq. © are displaced away from the track-ends, FinaUyj golving Eq . |A3] with boundary conditions Eq. \M 
as illustrated in Fig. |5J Eq. |6| then becomes 



we get 



(Al) 





I— > 





<$>P S 



i+1 



for < x < e, 
for e < x < I — e, 
for I — e < x < I. 



where e > 0. We recover the situation described in the (A5) 
main body of the text in the limit e 0. Eqs JA2] and [A5j together with the limit e 0, lead to 

In the periodic steady state we have Pf s (t+r) = Pf s (t) Eq. |19| 
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